HYDRODYNAMICS FOR A NON-CONSERVATIVE 
INTERACTING PARTICLE SYSTEM 



GLAUCO VALLE 



Abstract. We obtain the hydrodynamic limit of one-dimensional interacting 
particle systems describing the macroscopic evolution of the density of mass 
in infinite volume from the microscopic dynamics. The processes are weak 
pcrtubations of the symmetric exclusion process by shift operators describing 
the spread of particles around positions where new sites are created. 

1. Introduction 

The basic model of interacting particle systems that gives a rough microscopic 
description of the evolution of the mass density profile of an incompressible fluid is 
the exclusion process. Let p(-) be a finite range symmetric transition probability 
function on Z. The simple symmetric exclusion process (SSEP) on Z associated 
to p(-) is a Feller process with configuration space f2 = {0, 1} Z whose evolution 
can be described in the following way: initially each site of Z can be occupied or 
not by a particle, each particle of the system waits, independently of any other 
particle, an exponencial time of parameter one and at that time it choses another 
site according to p(-) and jumps to the chosen site if it is unoccupied. It is well- 
known that the SSEP has a hydrodynamic behavior under diffusive space-time 
scaling whose hydrodynamical equation is Laplace's equation, see jl]. 

Consider the following alternative description of the exclusion process: each site 
is occupied by two types of particles, black and gray particles as in figure [T] each 
black and gray pair of particles, independently of any other such pair, exchange 
positions after an exponencial time of parameter p(x — y) where x is the position of 
the black particle and y is the position of the gray one. In this way, black particles 
are associated to ones and gray particles to zeros for a configuration in f2. 
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Figure 1. A configuration for the exclusion process berore and 
after one extra black particle has entered 
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We describe now a pertubation of the exclusion process where new black particles 
are allowed to enter the system. The idea is to push a black particle between two 
nearest neighbor particles triggering an spread of particles centered in the position 
of the new particle, see figure [l] But after particles are translated, they are located 
in Z + | instead of Z. To deal with this behavior, we enlarge the configuration 
space to f2 = N x (r x U T 2 ), where 

ri = {r? 6 {0,l} zu ( z +s) : n(x) = forall i£Z+|) 

and 

r 2 = [v G {0, l} ZLJ ( z +5) : T}{x) = forall x G zj 

for Z+ g = {x+ g : .t G Z}. The system alternates between a SSEP on Z and a SSEP 
on Z + I , where the interchanges are a result of a superposed dynamics describing 
the spread mechanism that opens space for one new particle. This dynamics is 
described as follows: 

(i) Let h : [0, T] x K — > R + be a smooth bounded function such that there 
exist constants C > and /? > such that 

ft(t,«) < Ce"" 1 " 1 , (1.1) 

for every u G K and t > 0. In particular, 

sup C(t) < 00, (1.2) 
te[o,T] 

where 

/+00 
h{t,u)du, 
-00 

which is a bounded Lebesgue measurable function. 

(ii) Fix A^ > as the scaling parameter and let s denote a fixed macroscopic 
time. If the state of the system at macroscopic time s is (n, rj) G Ct with 
rj G (resp. rj G r 2 ), then we have that independently at each site 16Z 
(resp. x G Z + i) and at rate /i(s, x/AT) the system changes to (n + 1, ry') 
with ry' G r 2 (resp. Ti). 

(iii) The configuration ?/, mentioned in (ii), is obtained from r\ in the following 
way: all the system strictly at the right of site x (x included) is translated 
by 1/2 units to the right; all the system on its left hand side is translated 
by 1/2 units to the left; we set rf(x) = and rf(x— |) = 1. 

We are only going to consider initial cofigurations on f2 of the kind (1,7/) for 
rj G Ti, and therefore the system is always in 1^ for n odd and in T 2 for n even. 
Thus as mentioned before during each interval of time where n is od d (re sp. even) 
the system behaves as a SSEP on Z (resp. Z + |). From condition (1.2 1, we only 
have a finite increase of mass in finite time which is required for the system to 
be well defined. The system we obtain will be called the Exclusion Process with 
centered spreading mechanism (EPCS). 

In order to prove the hydrodynamic behavior of the EPCS we perform two 
transformations on the configuration space obtaining a system we can deal with 
more easily. The first transformation is obtained by reversing the roles between 
zeros and ones. The second transformation makes the process into a SSEP on 
Z with a superposed dynamics. To describe this dynamics set as the scaling 
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parameter and s as a fixed macroscopic time. Let b : [0, T] x K — > K + be related to 
h(-,-) by the equality 

b(t,u) = h(t,u- J o ^Y-d s y (1.3) 

Then for a site i 6 Zat rate 6(s, x/N), independently of any other site, the particles 
at the right hand side of x are translated one unit to the right. 

The aim of the first transformation is to obtain a system for which we do not 
need to deal with incoming particles. The second transformation is to simplify the 
configuration space obtaining a gradient modified exclusion process. The system 
we get will be called the exclusion process with right sided spreading mechanism 
(EPRS). 

The intuitive reasoning behind ( 1.3 ) is that the EPRS is obtained from the EPCS 
by translating the system one half units to the right each time mass increases by 
one unit and thus b at time t should be the translation of h in the space variable 
backwards by a half of the total mass that has entered the system before time t. 
Futhermore 

n+OO 

b{t,u)du = C(t), (1.4) 



and from (1.2 1 and the formal description of the EPRS through its generator we 
can show that the process is well-defined, see Section [3] 

The EPRS is a weak pertubation of the SSEP by shift operators and to prove its 
hydrodynamic behavior we have to deal with the particularities of this pertubation. 
The proof presented here is based on the entropy method of Guo, Papanicolaou and 
Varadhan [3] and an adaptation by Kipnis, Olla and Varadhan [3] to weak pertu- 
bations of the SSEP, but in [9] the pertubation of the Kawazaki dynamics is by 
Glauber instead of the shift operators for the EPRS. Hence this paper adds to the 
efforts to extend the results on hydrodynamics to new classes of non-conservative 
systems. Concerning the proof of the hydrodynamic behavior of non-conservative 
interacting particle systems, we have that similar arguments have also been em- 
ployed in [11] . Other papers as [6] and [7] also deal with adaptations of the entropy 
method and in [lOj the relative entropy method is applied. 

The EPCS and its hydrodynamic behavior are described in Section [2] and the 
EPRS and its hydrodynamic behavior in Section [3] The last sections, Sections [4] 
and [5] are devoted to the proofs of the hydrodynamic limits for the EPCS and the 
EPRS respectively. 



2. Hydrodynamics for the EPCS 

The EPCS is formally described through its generator which is an operator on 
the space of local real functions with domain Q,. Denote it by 



■>N 



N 2 L + L"' s 



The operator L is related to the SSEP dynamics associated to a finite range sym- 
metric transition probability p(-) and is equal to 

1 ver 1 ^2^2p(z){L x . x+z + L x+ZyX } + l ner2 ^2p(z){L x , x+z + L x+ZyX } 

xez Z £Z zgz+i-zez 
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where, for every local function F : f2 — * E we have 

L x<y F(n,ri) = r,(x)[l - V (y)} [F(n, r,^) - F(n,r,)], 

for every (71,77) G ^ an d x, y > 1, with rf' y being obtained from configuration 77 
with spins at x and y interchanged: 

)v(y) > i£z = x 
r)(x) , if z = y 
r)(z) , otherwise . 

The operator L^' s is related to the spread of the mass around the position where 
a raindrop falls: 

f N <» — 1 Vr^ s j-1 \ "* fN, s 
J-> r — i^er! 2^ W,x "+" - L ^er 2 ^ 

where, for every local function F : — » E 

77) = /i(s, a;/7V)[F(n + 1, f x r?) - F(n, 77)], 
where, for every y G ZU(Z+ |), 



(TxV)(y) 



v(y+l), if y < x- \ 

1 , if y = x - 2 
, if y = x . 



Note that condition (1.4) implies, as an application of the Borel-Cantclli Lemma 
that in finite time we only have a finite increase of mass. Thus, from the results 
presented in [8] there exists a non-homogeneous Feller process associated to the 
generator which we are going to denote by (rig , fj^ ) s >o- 

2.1. Weak solutions of the EPCS hydrodynamic equation. The Hydrody- 
namic equation for the EPCS is the convective diffusion equation 

d t p(t, u) = a 2 Ap(t, u) ~ §<9„(7(i, u)p{t, u)) + h(t, u), (t, u) G (0, T) x E, 
p(0, u) — po(u), u eR, 

(2.1) 

where 

/u r+oo y 

h(s,v)dv — / h(s,v)dv, a 2 = ~ z 2 p(z) . 
-00 Ju 2 zgZ 

In order to remove the external field h(t,u) from equation (2.1l we apply the 
transformation p = 1 — p. This is also required to obtain the hydrodynamical 
behavior of the EPCS from that of the EPRS, since from the microscopic point of 
view it is equivalent to reverse the roles between zeros and ones. Thus, we say that 
a bounded measurable function p : [0, T) x E — > E is said to be a weak solution of 
equation (2.1 1 if p : [0, T) x E — > E is a weak solution of 

d t p(t, u) = <r 2 Ap(t, u) - |5„( 7 (t, u)p(t, u)), (t, u) G (0, i) x E, 
p(0,u) = po(u), iieE. 1 j 

Denote by C^QO,! 7 ] x E) the space of real continuous functions with compact 
support on [0, T] x E which are continuously differentiable in the second variable and 
by C ' ([0,T] x E) the space of real functions with compact support on [0, T] x E 



HYDRODYNAMICS FOR A NON-CONSERVATIVE INTERACTING PARTICLE SYSTEM 5 



which are continuously differentiablc in the first variable and twice continuously 
differentiable in the second variable. Let po : R — > R be a fixed bounded measurable 
function. We say that a bounded measurable function p : [0, T) x R — > R is a weak 
solution of equation (2.2 1 if 

(a) p{t, u) is absolutely continuous in the space variable and d u p(t, u) is a locally 
square integrable function on (0, T) x R such that for all < t < T and for 
every function G G Cq 4 ([0,T] x R) we have that 



ds I duG(s,u)d u p(s,u) 



ds / dud u G(s,u)p(s,u) 



(b) For every function G G G 1,2 ([0,T] x R) and every i G [0,T] we have that 

duG(t,u)p(t,u) — / duG(0,u)p~o(u) = ds dud s G(s,u)p(s 7 u) 
Jw Jo Jr 

ds I —a 2 ( dud u G(s,u)d u p(s,u) + - I dud u G(s,u)j(s 7 u)p(s,u) 



Existence and regularity for equation (2.2), and thus for equation (2.1 1, follows 



from existence and regularity for equation (3.1 1, this is the contend of Lemma 4.2 
The required uniqueness result is given by theorem 5.2 in [5]. 



2.2. The EPCS hydrodynamic behavior. We denote by "P(fi) the set of proba- 
bility measures on f2. Let v a G Vifl) be the Bernoulli product measure of parameter 
a G [0, 1] on f2. Given any two probabilities p, v in V(Q), we denote by H{p\v) the 
relative entropy of p with respect to v. 



H(p\v) = sup jy /d/x - log J ( 



:^d^ 



where the supremum is carried over all bounded continuous real functions. 

A sequence of measures (pn)n>i in "PiP) ' IS associated to a initial profile po : 
M — ► M if for every 8 > and every continuous function G : R — * M with compact 
support 



lim a 



i^G(^/iV)^)- f duG(u) Po (u) 



> S 



0. 



Let D(R_|_, 17) denote the space of right continuous functions with left limits on 
fl endowed with the Skorohod topology. For each probability measure m on O and 
N > 1, denote by the probability measure on Cl) induced by the Markov 

process (n^ , r]^) s >o with generator and with initial measure m. 

For a measure p in V(£l), let T{p) be the measure on f2 induced by p through 
the transformation £ G fi — > (1, J?) G fi where 77 is the configuration in Ti which is 
equal to £ on Z. 

The hydrodynamical behavior of the EPCS is given by the following result: 

Theorem 2.1. Fix a sequence of initial probability measures (pn)n>i inV(£l) as- 
sociated to an initial profile po : R — * R+, bounded above by 1, such that H{p N \v a ) < 
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CN for some a £ (0, 1). Then, for any continuous function G : 
pact support, any 8 > and < t < T 



with 



lim , 

jV-«x> 1 W I 



i £ G{z/N)th{z) - [ duG(u) P (t, 



zezu(z+i) 



> <5 



where p is the unique weak solution of the convective diffusion equation (2.1) 



3. Hydrodynamics for the EPRS 

The EPRS is formally described through its generator which is an operator on 
the space of local real functions with domain Q. Denote it by 



N L + L 



N.s 



Here L is related to the motion of particles in a SSEP associated to a finite range 
symmetric transition probability p(-): 



L — ~y^^p(z){L x x+z + L x+zx } 



where, for every local function F : — > R and every x,y > 1, 

L XtV F(o = s(x)[i - av)} \ne' v ) - 

and £ x,y is the configuration with spins at x and y interchanged: 

!t(y) > if z = x 
^(x) ,Xz = y 
£(z) , otherwise . 

The operator L^' s is related to the translations of the system at the right of a given 
site: 



L 



N,s 



N,s 



where, for every local function F : CI — > R 
with 

{£(a; — 1) , if x > z , 
, ifs = z, 
, if a; < z . 

We refer to Liggett 's book [8], to a proof that £^ is the generator of a Feller 
semigroup. Therefore, there exists a non-homogeneous Feller Process associated to 
the generator £^ which we are going to denote by (^) s >o- 
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3.1. Weak solutions of the EPRS hydrodynamic equation. The hydro-dy- 
namic equation of the EPRS is also a convective diffusion equation: 

d t ((t, u) = <J 2 A((t, u) - d u (a(t, u)C(t, u)), (t, «)el + xl, 
C(Q,ti)=Co(«), ueR, 

where 

/x y 
b(s, u)du , a 2 = - z 2 p(z) . 
-°° zez, 

For a fixed bounded mesurable function Co : R — ► R, a bounded measurable 
function £ : [0, T) x R — ► K is said to be a weak solution of equation (3.1 ) if 

(a) £(t, it) is absolutely continuous in the space variable and d u ((t, u) is a locally 
square integrable function on (0, T) xK such that for all < t < T and for 
every function G G C^QC^T] x M) we have that 

ds I duG(s,u)d u ((s,u) = — / ds dud u G(s,u)((s,u) . 



(b) For every function G G C 1,2 ([0,T] x R) and every t £ [0,T], 

duG(t,u)^(t,u) — / G(0, tt)Co('u) = ds dud s G(s,u)(^(s,u) 
Jr Jo Js, 

ds < — a 1 \ dud u G(s,u)d u (^(s,u) + / dud u G(s,u)a(s,u)^(s,u) 



Existence and regularity for equation (3.1 1 follows from the proof of the hydro- 



dynamical behavior of the system described in Theorem |3.1| on Section |3.2| The 
required uniqueness result is given by theorem 5.2 in [5]. 



3.2. The EPRS hydrodynamic behavior. Let £>(R + ,f2) denote the space of 
right continuous functions with left limits on f2 endowed with the Skorohod topol- 
ogy. For each probability measure pL on and N > 1, denote by P^ the probability 
measure on £)(R_|_,fi) induced by the Markov process £j with generator £^ and 
with initial measure fi. Recall also the definitions of v a , the relative entropy from 
Section |2.2| and of a family of probability measures on Q, associated to a density 
profile. 

The hydrodynamical behavior of the EPRS is given by the following result: 

Theorem 3.1. Fix a sequence of initial probability measures {h n )n>i on il associ- 
ated to a initial profile Co : R — ► R, bounded above by 1, such that H(fj, N \is a ) < CN 
for some a G (0, 1). Then, for any continuous function G : R — > R with compact 
support, any 8 > and < t < T 



lim P^v 



Lj2G(z/N)Z t (z) - f duG(u)t(t,: 



> 6 







where £ is the unique weak solution of (3.1). 
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4. From the Hydrodynamic for the EPRS to the hydrodynamic for 

THE EPCS 



This Section is devoted to the proof of Theorem |2.1| The proof is based on 
a coupling and thus comparative arguments between the EPCS and the EPRS, 
therefore it follows from Theorem 3.1 which is shown in Section [5] For the coupling 
we use an auxiliary process which is obtained by translating to the left the EPRS 
by a half times the number of translations of the system that have already ocurred 
and then by reversing the roles between zeros and ones. In this way, if we put in 
correspondence the number of particles entering the EPCS with the translations 
in the EPRS, we have that the auxiliary process is basicaly the EPCS except for 
corrections in the positions of particles that are negligigle in probability as N — > 
+oo, where N is the scaling parameter. 

Recall from the statement of Theorem |2 . 1 1 that we will be considering a family of 
probability measures (/j, n )n>i on satisfying some conditions imposed there. Also 
from Section |2.2[ this family can be regarded as a family of probability measures 
on O through the correspondence between T± and f2. Thus it represents the initial 
distribuition for both the EPRS and the EPCS in this section. 

Aiming at the strategy described above, let us first make some considerations 
about the number of translations of the EPRS. At scaling N and for each x £ Z, we 
denote by (W^' N : < t < T) the independent non-homogeneous Poisson processes 
with rate b(s,x/N), < s < T, which represents the number of translations to the 
right of the EPRS ocurring from site x before time t. Next lemma relates in a weak 
sense the total number of translations of the EPRS with the macroscopic rate at 
which the system is translated. 



Lemma 4.1. For every S > and < t < T we have that 



lim 

N— >oo 



jlV 



N 







f C(s)ds 


> s 


10 





Proof: We show that for every u £ 



N 



lim P„ 



x< [uN] 



W; 



x,N 



h(s, v) dvds 



> S 



= 0. 



It is clear in the proof that we may replace u by +oo in the previous expression to 
get the limit in the statement. Since convergence in probability to a constant fol- 
lows from convergence in distribution, we only have to show that J2x<[uN] W*' N 

converges in distribution to L /" h(s, v) dvds. We are going to estalish the con- 
vergence of the appropriate characteristic functions. Since, for a Poisson process 
W with time varying rate A(s) we have that 



E[eMnW t }} = exp {(e^ - 1) J \(s)ds}, 
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we get from independence that 



E 



E »'■""} 



x<[uN] 



Tlx< [uN]E 



exp 



n x <[„^]exp|(e^ - 1) j h(s,x/N)ds} 
expj Y (e^ - 1) y /i(s,a;/JV)ds|, 



x<[uN] 



Now 



for every 7 G 

V (e&-l) / h(s,x/N)ds -> 17 / / 

w An JO JO J-i 



/i(s, v)dvds , as ./V — > 00 . 



K<[«iV] 

and then 



exp < 17 



{< 



/i(s, it)duds| , as iV 



"Ms E 

a;<[uAT] 

for every 7 > 0. Therefore we have proved the required convergence. □ 

Let (£,^)s>o denote the EPRS. To obtain the auxiliary system we consider the 
transformation from _D(R + ,f2) to D(M.+ ,Ci) which associates to the configuration 
^ the configuration (n^,f)^) given by 



and 



Thus 



1 - 




t»(x + *f), if (* + 



G Z 



otherwise. 



i ^ G(x/N)^(x) 



is equal to 



*E^( 



a;eZU(Z+i) 



- JV 



- JV 



N 



1 / 



-2 



if/ 



(4-i) 



By Theorem |3.1| and Lemma |4.1| the term in the rightmost side of the previous 
equation converges in probability to 

C(r) 



G [u 



-dr \ (1 — C( s i u))du 



(4.2) 



is the weak solution of (3.11. We 



for all < s < T, where ( : [0,T)xR- 
want to show that the same holds for the EPCS replacing the auxiliary process. 
To show this assertion we couple the EPCS with the auxiliary process in a proper 
way. Our coupling will allow first and second class particles on both systems. The 
difference between them is that first class particles have priority over second class 
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particles this means that if a first class particle attempts to jump over a second 
class particle then they exchange positions. On the other hand, if a second class 
particle attempts to jump over a first class particle the jump is suppressed. The 
description of the coupling follows below: 

(i) Both systems start with the same configuration. For both systems, particles 
in the initial configuration are labeled in an ordered way from left to right. 
All these particles are first class, indeed only the first class particles are 
going to be labeled. 

(ii) A jump attempt of a single first class particle occurs for the EPCS if and 
only if a jump attempt in the same direction occurs for the first class particle 
with the same label in the auxiliary process. Second class particles on both 
systems move independently of any other particle. 

(iii) The appearance of a particle and the induced translation may occur in 
only one of the systems at a time or in both of them simultaneously Recall 
that for the EPCS a particle appears at site x — \ at time s with rate 
h(s,x/N) — b(s, + J Q S ^p-dr) and for the auxiliary process with rate 

H s ) J? + 2J\r)- I n order to deal with the fact that a coordinate process 
might be evolving on Z while the other is evolving on (Z+ \ ), define x' = x 
if both coordinate processes have received the same amount of new particles 
and are evolving on the same lattice, and x' = x + ^ otherwise. We have 
one of three possibilities: 

(1) If h(s, x/N) > b(s, + only the EPCS is translated around site 

' n N 

x at time s at random rate h(s,x/N) — b(s, ^ + ^) with a particle 
appearing at site x — \ . This particle becomes a second class particle 
for the EPCS. 

(2) If h(s,x/N) < b(s, + 2jf), only the auxiliary process is translated 
around site x' at time s with a particle appearing at site x' — \ at 

random rate b(s, ^+2n) — M s > x/N). This particle becomes a second 
class particle for the auxiliary process. 

(3) Both systems are translated simultaneously at time s, the EPCS around 
site x with a particle appearing at site x— \ and the EPRS around site 
x' with a particle appearing at site x' — ^ . This happens at random rate 

' n N 

min{ft,(s, x/N),b(s, ^ + jn)}- These particles are first class. Since we 
have new first class particles, we have to relabel them in a way to keep 
the particles ordered from left to right. Let us explain how we relabel 
these particles: before the transition, due to (1) and (2), we may have 
pairs of equally labeled particles that do not share the same positions. 
Let the positions of one of these pairs be denoted by (z, w) with the 
first particle belonging to the EPCS and the other belonging to the 
auxiliary process. If z < x and w < x' or z > x and w > x' then the 
two particles remain to share the same label after transition. Now sup- 
pose that there exists a number I (necessarily finite) of such pairs with 
positions (zj,Wj), 1 < j < I, such that Zi < Zj < x and x' < Wi < Wj 
or x < Zi < Zj and Wi < Wj < x 1 , for every i < j. Now we consider 
new pairs (zi - \,x' - \), (z 2 - \,w\ + §), ... , (zi - \,wi-i + \), 
(x — |, W[ + |) in the case Zi < Zj < x and x' < Wi < Wj or new pairs 
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Or- 



|), (zi + i 



w 2 



2), ... , (Zl-1 + 2> W l- o)' 



I) 



in the case a; < z, < z„ and w, 



< Wj < x' . 



Finally give the same label 



to both particles in each of these pairs of positions. Note that in this 
transition the distance between two equally labeled particles can only 
decrease. 

Following this we obtain a Markov process on f2 x which will be represented by 
((71^,77^), (n^ , fj^)) s >o where (n^ ,r)^) is distributed as the EPCS and (nf ,fj^) 
as the auxiliary process, this will be called the coupled system. For a formal de- 
scription, one can write explicitly the generator of the coupled system or use Harris 
graphical construction which are already standard in the literature, we refer here 
to Liggett 's book [8]. 

Note that the maximum distance at time s between two equally labeled particles 
in each one of the two coordinate processes of the coupled system is bounded by 
the total number of second class particles. Let this random number be denoted by 
Jf . We have that (J^) s >i is pure birth process with random rate 



X 

N 



C{r) 



dr 



bis, 



N 



2N 



By Lemma 4.1 and the cxponcncial bound on b given in ( 1.1 1, the integral in the time 



interval [0, T] of the previous expression divided by N goes to zero in probability 
as N — > oo. Therefore, for every e > and < 7 < 1 we can fix Ng = No(e) 

J N 

such that for N > No the random variable -jj- is stochastically bounded on a set 
of probability greater than 7 by a Poisson distribution with mean eN and scaling 
parameter iV _1 . This Poisson distribuition converges to zero in probability since 
it has characteristic function exp{eiV(e^ — 1)} which goes to 1 as N —> 00 and 

jN 

0. Hence also converges to zero in probability as N — > 00. Note that 



n 
N 



N 



- n 

DO 



N] 



< , and therefore 



N 



also converges to zero in probability as 



Now as in the description of the coupling we label the first class particles of 
the EPCS in a ordered way from left to right by denoting their positions at time 
s by (s), for k G Z, and we do the same for the auxiliary process denoting 
the positions by XP(s). We do the same with the second class particles denoting 
their positions respectively Yjr (s), 1 < k < m, and Yjf(s), 1 < k < I, where 
m = m(N, s) is the number of second class particles in the EPCS and I — l(N, s) 
is the number of second class particles in the auxiliary process. In particular, 
Jf = m + l. 

Finally we are in condition to compare integrals with respect to the empirical 
measures associated to 77 and t). We have that 



1 ]T G{x/N)r£{x)-± J2 G(x/N)^(x) 



x6ZU(Z+J) 



(4.3) 



is equal to 



fcez 



N 



G 



N 



1 

N 



k=l 



fe=l 



TO - 
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Since G is a smooth function with compact support we have that this last expression 
is bounded above by 

c (G )(^ ^»-^" + f), 2C(G) f, 



where C(G) is some constant depending on G. Therefore (4.3 1 converges to zero in 
probability as N — > +00 and from (|4.1|) and (|4.2h any (5 > and < < < T 



lim , 



i J] G(z/N) m (z) - I duG(u)p(t,u) 



zezu(z+i) 



> <5 



= 



where p is given by 

u) = 1-C(t,u + J ^^ ds ) > for ever y (*> M ) G [0, T] x M . (4.4) 



Now with the next lemma we end the proof: 



Lemma 4.2. Let ( : [0, T) x K - 
1 — p(-, •), wii/i j o(-, •) defined in (4-4K * s a weafc solution of (2.2) 



be a weak solution of {3.1). Then p(-,-) 



Proof: Note first that for a function G G Cq' ([0,T] x E), the function F defined 
by 

F(t, u) = G(t,u- J ^y^rfs) . for every (t, u) G [0, T] xR, (4.5) 

is also in Cg' 1 ([0, TjxR). Thus we have that (a) in the definition of weak solutions 
of equation (|2.2| follows for p directly from (4.4) and (a) in definition of weak 
solutions of ( |3.1| by a change of variables on both sides of the integral. 

Therefore it remains to verify that (b) in the definition of weak solutions of 



equation (2.2 1 follows from (b) in the definition of weak solutions of equation (3.1 1 
Fix a function G G C ^ 2 ([0,T] X R), we have that 



duG(t,u)p(t,u)- / duG(0,u)p o (u) = / du F(t,u)((t,u) - / du F(0, u)Co(u) , 



(4.6) 

where F is as in (4.51. We can apply the equality (b) in the definition of weak 
solutions of (3.1l if F G C ' ([0, T] x R), however we cannot guarantee that the 
first derivative in the first variable is differentiable since 



d t F(t,u) = d t G(t,u 



C(s) 



ds) - 



C(s) 



t.u — 



C(s) 



ds 



and C is not necessarily continuous. To overcome this problem, we rely on some ap- 
proximation results. Let D : [0,T] — > M. + be the function D(t) = | J Q * C(s)ds. The 
function D is a continuous increasing function with bounded, thus integrable, deriv- 
ative. By Theorem 1 in Section 4.2 of [2], there exists a family of smooth functions 
D £ , for e > 0, obtained from D by the convolution with standard mollifiers, such 
that D e converges uniformly to D and dD e converges to dD in L 1 ([0, T]). In accor- 
dance to (4.5 ), let F e be the smooth function F e (t, u) = G(t, u—D e ), for every (t, u) G 
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[0,T] 



Since F e € C 1,2 ([0,T] x K) then 



duF e (t,u)C(t,u)- / duF e (0,u)( o (u) = / ds / dud s F e (s, u)C(s, u) 



ds I -a' 



du d u F' i (s, u)d u ((s, u) + / dud u F' i (s,u)a(s,u)((s,u) 



Now, since ^ is bounded, <9„£ is locally square integrable and from the definition 
of F e , we can stablish by straightforward manipulations the convergence term to 
term as e — > in the right hand side of the previous equality. Therefore 



duF(t,u)((t,u) - / duF(0,u)( o (u) = ds dud s F(s,u)((s,u) 



ds < -o' 



dud u F(s,u)d u ((s,u) + / dud u F(s, u)a(s, u)((s, u) > , 



where the right hand side can be rewritten as 



ds / dud s G(s,u)p(s,u) + / ds < — a / dud u G(s,u)d u p(s,u) + 



dud u G(s,u) a(^s,i 



*C(a) 



ds 



C(s) 



p(s,u) 



(4.7) 



Moreover we have by a simple computation that 



a[ s,u + 



C(s) 



ds 



C(s) _ 7 (t,u) 



and thus from (4.6 1 and (4.7 1 we have that (b) holds in the definition of weak 
solutions of equation (2.2). □ 



5. The proof of the Hydrodynamic limit of the EPRS 



Denote by M. = A4(R), the space of positive Radon measures on R endowed 
with the vague topology. Integration of a function G with respect to a measure 
7r in M. will be denoted (ir,G). To each configuration £ £ f2 and each N > 1 we 
associate the empirical measure ir N — 7r Ar (^) in VW, where 

zGZ 

Let D([0,T], M.) denote the space of right continuous functions with left limits on 
M. endowed with the Skorohod topology. For each probability measure /j, on f2, 
denote by the probability measure on D([0,T},M) induced by and the 
empirical m easu re ir N 
Theorem 



3.1 



states that the sequence Q N converges weakly, as N — > oo, to the 
probability measure concentrated on the absolutely continuous trajectory ir(t, du) = 
((t,u)du whose density is the unique weak solution of (3.1l (see [1]). The proof 
consists of showing tightness of Q N N , that all of its limit points are concentrated 



on absolutely continuous paths which are weak solutions of (3.1 ) and uniqueness of 
solutions of this equation. 



We have already discussed uniqueness of weak solutions of (3.1 ) in Section 3.1 
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Note that all limit points of the sequence are concentrated on absolutely 
continuous measures since the total mass on compact intervals of the empirical 



measure 7r is bounded by the size of the interval plus 1/N. 



In order to show that all limit points of the sequence 



weak solutions of (3.1 1 we will need the following result: 



are concentrated on 



Lemma 5.1. For every smooth function G : [0,T] x 
and S > 

ft 



with compact support 



lim sup i 



t,7V 



N- 



sup 

0<t<T 



5.1 



«,G)-«,G)- / {(n»,d s G)+ 
Jo 

+ a 2 <7rf , AG) + (Trf ,a s VG)}ds\ > 5] = . 



ft TV 



We shall divide the proof of Theorem |3.1| in four parts: We start proving tight- 
ness in Sect ion |5.1| The Section |5.2| is devoted to prove of Lemma |5.1| From 
Lemma 



to conclude the proof that all limit points of the sequence 



concentrated on weak solutions of (3.1), we have to justify an integration by parts 



to obtain conditions (a) and (b) in the definition of weak solutions of (3.1 1. This is 
a consequence of an energy estimate which is the content of Section 5.3 



5.1. Tightness. The sequence is tight in the space of probability measures 
on -D([0, T], A4), if for each smooth function with compact support G : R — > R, 
(tt^,G) is tight as a random sequence on £)([0,T],R). Now fix such a function, 
denote by Tt — c 2 (-7r s : s < t), t > 0, the natural filtration on D([0, T], Ai), and by 
Tt the family of (^-"t)-stopping times bounded by T. According to Aldous [I], see 



also Chapter 4 in [4], to prove tightness for 
two conditions: 



, G) we have to verify the following 



(i) The finite dimensional distributions of 

(ii) for every e > 



, G) are tight; 



lim limsup sup sup [|(7rf ,G) - (tt? +9 , G)\ > e] = 



(5.1) 



N^oo tET t 6K7 

Condition (i) is a trivial consequence of the fact that the empirical measure 
has finite total mass on any compact interval. In order to prove condition (ii), 
Let us first consider for each smooth function H : [0, T] x M — > R the associated 
(J-t)-martingale vanishing at the origin 



M. 



H.N 



H) 



H) 



{d s + C»){K»,H)ds 



(5.2) 



,H.N 



In (ii) the function G does not depend on the time, however the martingale M t 
is defined for functions varying on time, because we will need it later in the proof 
of lemma 5.1 An elementary computation shows that (d s + £^)(tt n , H) is given 



explicitly by 



(n N ,d s H) + (n N ,A p N H) + ^ £ U £ b(s, x/N) \ X7 N H(s, z/N)^z) , (5.3) 
where A^ and Vjv denote respectively the discrete Laplacian and gradient: 



A p N H(s, x/N) = — £p(«){ff(«, (a; + z)/N) + H(s, (x - z)/N) - 2H(s, x/N)} , 
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V N H(s,x/N) = N{H(s,(x+l)/N) - H(s,x/N)}. 
The quadratic variation of M t ' is denoted by {M H ' N ) t and explicitly given by 

f x V y — 1 ^ 

/ <Ma^5> (z) E ( e v W ff( SlW /jv))Vwii-e(!/)]+ 

z6Z + \x-y\=z u>=xAy 

+ ^ E a ( s ' z/N)VnH(s, z/N)*g(z) 

Z 

l - ^ A w/N)V N H(a, z/N)V N H(s, w/N)g(z)g(w) (5.4) 



TV 3 



^ avw — J- 

"^E{ E 6(*,a:/JV3}ff(a^Vti;/JV)V^( S> zA«;/JV)^(2)^(«;)}, 

where a N (s, z) = N^ 1 J2 x <z X /N) . Therefore, for a G not depending on time 
as before 



G) « , G) = - + I £f « , G)ds 



From the previous expression and Chebyshev inequality, (ii) follows from 

lim limsup sup sup E* N \ M°r% - M?' N } = (5.5) 

A^oo t£T t 9<7 L J 

and 



lim limsup sup sup E^jv 

7^° N-nx> tET t 9<7 



r+e 



= 0. 



(5.6) 



Now we show (5.5 1 and (5.6 1, completing the proof of tightness and finishing this 
section. 



Proof of (5.5 ): From the optional stopping theorem and the martingale property 
\(M^ - M?- N Y] = E£ [(M G ' N ) T+e (M G ' N ) T ] . 



E 



N 



Hence, applying formula (|5.4|, by the Taylor expansion for G, we have that 



E 



N 



{M GN_ M c iNf 



< 



C{G,b) 
N 



(5.7) 



and (|5.5|) holds. Indeed, to deal with the last term in (|5.4|), use the fact that b is 

lal: 

2 



Ricmann integrable together with the inequality 
1 ^ zVw-l _ 

y: E E E kw*) 



iV 3 ^ ^ iV 3 

\z\,\w\<CN X = ZAW + 1 \z\,\w\<CN W<X<Z 

z^u, z>n, 

9 1 CN 

^ E {ir E «'.*/")}■ 



□ 



Proof of (5.6): From formula (5.3 1 and the Taylor expansion for G we obtain that 

i 

£f(7rf, G)ds <C(G,6)0. □ (5.8) 



E 



N 



T + 9 
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5.2. The proof of Lemma ( pTT] ). Let G : [0, T] x K -> M be a smooth function 
with compact support. By Doob inequality, for every S > 0, 



sup |Mf ,iV | > <5 

0<t<T 



< 4S~ Z E 



2tu,N 



{M^ N f\ = 4<T 2 E^ [{M G > N ) T } , 



G iV 

which, by the explicit formula for the quadratic variation of M t ' , see the proof 
of ((5T5}, is bounded by C{G, b) ON' 1 . Thus, for every S > 0, 

= 0. 



lim P^v 



sup |Aff * w | > S 

0<t<T 



(5.9) 



Using ( 5.3 1 to expand the martingale expression in ( 5.2 ) . Since the Taylor expansion 

C(G) 



gives us that 

sup \N[G(s,x + l/N) - G(s,x/N)} - VG(s,x/N)\ < 



se[o,r] 



N 



and A^H(s, x/N) is equal to 

1 \ -> , v r t x — y + w 



x — y + iu — I 
TV 



)>■ 



we may replace and V^v in ( |5.9| by the usual laplacian and gradient. More- 
over the exponencial bound on b given at ( |1.1[ ) implies that the Riemann sum 
iV -1 J2 x <[uN] b( s i x /N) converges as N — > 00 to J™ b(s,v)dv — a(s,v) and we can 
also replace one by the other in (5.9). Thus we obtain that 



lim 



sup 

0<t<T 



(ir? ,G) - (ir$ ,G) - / ds{(^,d s G)+ 
Jo 

+a 2 (7rf,AG) + (7rf,aVG)}| > 6} =0 



for all 5 > 0. □ 

5.3. An energy estimate. The next result allows an integration by parts in the 
expression inside the probability in the statement of Lemma 5.1 proving Theorem 

i3~n 

Proposition 5.2. Every limit point of the sequence Q N is concentrated on paths 
£(t,u)du with the property that C(fi u ) * s absolutely continuous whose derivative 
d u ((t,u) is in L 2 ([0,T] x R). Moreover 



ds / du H(s ) u)d u C,{s,u) = — I ds / dud u H(s 1 u)(^(s,u) , 
Jr Jo Jm 

for all smooth functions U:[0,T]xR^R with compact support. 

To prove the previous theorem we make use of the following energy estimate: 

Lemma 5.3. There exists K > such that if Q* is a limit point of the sequence 
Q N then 

( r T r r T r 1 

<K, 









Eq, 


sup < 








JO JTS 



ds j dud u H(s,u)((s,u) — 2 / ds duH(s,u) 2 ((s,u) 

Jo Js, 

where the supremum is taken over all functions H in Cq :1 ([0,T] x R). 
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Proof of Proposition 5.2 Let Q* be a limit point of the sequence Q N . By 
Lemma 5.3 for Q* almost every path ((t,u) there exists B = B(Q > such that 

[ ds [ dud u H(s,u)((s,u) -2 [ ds [ duH(s, u) 2 < B, (5.10) 
Jo Jr Jo Jr 

for every H S O ' ([0, T],R). Note that, since £ < 1 we were able to suppress it in 
the last integrand. Equation (5.10) implies that 



X(H):= [ ds [ dud u H(s,u)((s,u) 
Jo Js. 



is a bounded linear functional on Cq ,:l ([0,T],E) for the L 2 -norm. Extending it to a 
bounded linear functional on L 2 ([0, T], E), by Riesz Representation Theorem, there 
exists a L 2 function, denoted by d u ((s, u), such that 

X(H) 



ds / du H (s, u)d u ((s, u) , 
Jr 

for every smooth function H : [0,T] x E — * K with compact support. □ 

To show Lemma [5^3] consider the space C° A ([0,T] x E) endowed with the norm 

00 

||ff||o,i = ^2 n {\\H 1)}IU + ||a u fl"l{(n,n + l)}|| 0o }. 



n=0 



Proof of Lemma |5.3| For every e>0, <5>0, i?:E^E smooth function with 
compact support and £ € {0, 1} Z , denote by Wn(c, S, H, £) the following expression 

00 00 [eN] 

J2 H( X /N) — {t; SN (x)-t SN (x+[eN])}~- ]T H(x/N) 2 — ^ N {x + y), 



y=0 



where £ s (x) = S 1 Yly=x £(?/)■ We claim that there exists K > such that for any 
dense subset {ff, : I > 1} of C^QO, T] x E) , 

dsW N (e,8,H i (s,-),Z s )'> <K. 



N 



lim lim E,, 

5^0 iV->oo M 







max < 


Jo 


l<i<A: 





(5.11) 



for every k > 1 and every e > 0. We postpone the proof of (5.111 and we use 
it in the sequence. Since Q* is a weak limit point of the sequence Qn, which is 
concentrated on absolutely continuous trajectories, it follows that 

ds I du 



lim sup Eiq 
s^o 



max 

Ki<k 



(i 



u+S 



Csdv — S 



u+t+S 



Li + C 



^e^H^s.u) 2 / dvU- 1 < 

J u \ J V 



Cs dv 



Csdv' 



(5.12) 



< K. 



for every k > 1. Since, 



1 / duH(u) IS' 1 / (sdv-6- 1 



u+S 



u+e+5 



Cs dv 
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is equal to 



du 



H(u) - H(u - e) 



Cs dv } , 



letting S — > and then e — >• 0, it follows from (5.12) that 

rT 



Em* 



max 

Ki<k 



ds / du{d u Hi(s,u)C(s,u) - 2Hi(s,u) 2 C(s,u)} 



< K . 



To conclude the proof we apply the monotone convergence theorem for k — > oo to 
replace the maximum over 1 < i < k by the maximum over i G N in the integrand 
above, then note that 



ds J du {d u H(s,u)((s,u) — 2H(s,u) 2 ((s,u)} , 
is continuous as a real function on (Cq ,:l ([0, T] X K.), || • ||o,i)- D 



Proof of (5.11 ): Since H is a continuous function, an integration by parts justify 
the replacement of Wn(c, 5, H,£) as S — > in (5.11 1 by 



AT] 



^ 1 9 _ , 1 L J 



x— — oc> 



x= — oo 



which we denote by Wjv^e, ii", £). By the entropy inequality 



(5.13) 



E 

is bounded by 

B(ji N \u a ) , 1 



iV 







max < 


Jo 


l<i<k 









dsW N (e,Hi(s,-),Z s ) 



N 



N 



N 



exp < N max / dsWjv(e, i?i(s, •), 6) 

l<i<kj 



The bound on the entropy given in the statement of Thcorcm |3 . 1 | and the elementary 
inequality e max ai < ^ e ai imply that this last expression is bounded by 



C+-logE 



N 



Y, ex p\ N J o dsW N (e,H t (s,-),i s ) 



Here, since maxjlim sup w N 1 log ajv , lim supjy N 1 log is greater or equal to 
lim sup ^ N^ 1 logja^r + &at}, the second term is dominated by 



max lim sup — log E^ 

1<-Kfc N ^ N 6 "« 



cxpiN dsW N (e,Hi(s,-),Z s ) 



(5.14) 



We have that the previous expression is bounded by 
max f dssup { f W N (e, H t (s, •), £) - iV£>(/)| + j- £ 



1 ^ b(s,x/N) 
2(1-") 



(5.15) 
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where the suprcmum is taken over all cylindrical densities / with respect to v a , and 



D(f) 



fLy/ fdv a 



\ E p(*) / ^)[i-^+^)][\/I^)-vW)f^ 



is the Dirichlet form evaluated on /. We assume (5.151 and show it just after the 
proof of (5.11 1. We just have to estimate the first term in (5.151 whose limsup as 
N — > oo we are now going to show that is in fact non-positive. Write 

\t(x)-ax+[eN])}f(i)u a (d£) = 

2>1 j=0 J 

+ 2 E P(z) I {£(x + [z^eNlz) - £(z + [eN])}f(^ a (dO 



where the first term in the right hand side can be rewritten as 

[z -1 eAT]-l 



X+jz,X+(j-\-l)Z\ 



z >i j=a 
Hence, we have that 



oo „ 

E H(s,x/N) j{£(x)-ax + [eN])}f(Z)v a (dO 



is, for each B > 0, bounded above by 

E bJ2p(z) E 3(-*>i./) 



X= — OO 2I>1 

oo 



E 



H(s,x/N)' 



X — — OQ 



B 



[z~ 1 eN]-l 

E Z(+,x,j,f)+0(l), (5.16) 

2 >1 j=0 



where n(±,x,j,f) denotes the integral 

& + u + i)z)[i - ^ + jz)] | yjf^+j^+u+i)*) ± I ^(do- 



Multiplying (5.161 by (eiV) 1 , we bound it above by 

2 oo [.A] 

+ S E H(s,x/N) 2 — E I Z(x + y)f(0vc,(dt) + O((eN)- 1 ). 
x——oo y—0 

Now taking B = N we obtain from ( 5.13[ ) that 

W N (e, H(s, ■),0f(0v*(<%) < ND(f) + Oi(eN)- 1 ) . 
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Thus the limsup as N — > oo of the first term in (5.15) is non-positive and (5.141 is 
bounded by 

1 v-^ b(s,x/N) 

IV x£Z V ' 

which is finite since b is integrable. □ 

Proof of (|5.15[): Let U e N (s,£, s ) := W N (e,Hi(s, •),£„). We wish to estimate 



1 

N 



logE 



cxpj^ NU e N ( S ^ s )dsj 



(5.17) 



Define 



/(&_,) exp / NU»(s + r,Z r )dr 



for every bounded real function / on {0, 1} Z . We have that 



E 



A' 



exp<j / NU?(s,£ a )ds 



= I P Q N t l dv a <ll (P»lf dv, 



(5.18) 



We show below that 

D. f (-Ff.l) 2 < j 2s»P { / WDf(s, - JV 2 D(/)} + 



From (5.181, this inequality and the Gronwall inequality, we have (5.151. 



To complete the proof we show (5.191. By Feynman-Kac formula, {P s .t '■ < 
s < t} is a semigroup of operators associated to the non-homogeneous generator 
N 2 C^ + NU^ (s, •). Moreover, the first Chapman-Kolmogorov equation holds: 
d s P s , t = -(N 2 C^ + NU e N (s, -))P Stt . Hence 

1 



J s J (P» t l) 2 dv a = J (N 2 £» + NU t N (s, -))« t l) < 4 1 dv a 

= J N 2 L^(P»l)P s N t l dv a + J{N 2 L + NU»( S ,-)}(P s N t l)P s N t l 

We shall estimate separately each term in this expression. 
Claim 1: 

/ L?"{I%1) P»l du a < J2 h ^0^ J (^l) 2 dv a • 



Proof of Claim 1: Denote P* t l by h. We are considering an integral of the form 



E f Ks, x/N)[h(r x t) - h(0] MO dv a , 

x6Z 



where t x are defined in Section [3] Add and subtract to this the term 



(l-£(x))h 2 dv a 
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We end up with two terms, the first is 

f[h(r x 0-h(0] 2 dv a - \ ]T b(s, x/N) f h 2 dv a 

which is negative and may be neglected, and the other is 



Claim 2: 



J {N 2 L + NU e N (s, 0}(i$l) P s N t l dv a < 



< sup|y NU?(a,Z)ffiu a (dZ) -N 2 D(f) j J (P s N t l) 2 dv a . 

Proof of Claim 2: We have the bound 

J {N 2 L + NU?(s, -)}(P s N t l) P s N tl dv a < Tf J (P s N t l) 2 dv a , 

where is the greatest eigenvalue of the generator N 2 L + NU^ (s, •). By the 
variational formula (see Appendix 3 in [5]) is equal to 



sup^ / VU?{8,Z)f($v a {dt)-N*D[f) 



□ 



(5.20) 



From the claims 1 and 2 we have (5.19). □ 
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